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We consider the Anderson localization problem on the infinite regular tree. Within the localized 
phase, we derive rigorous upper and lower bounds on the free energy function recently introduced by 
Aizenman and Warzel. With an additional but yet unproved regularity condition on this function, 
this leads to bounds on the critical disorder that induces localization. The latter are particularly 
useful in the limit of large connectivity where they match and confirm the early predictions of 
Abou-Chacra, Anderson and Thouless. 
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I. INTRODUCTION 

The strong effect of disorder on the transport properties of a quantum particle was first brought 
to light by the seminal work of Anderson Since then, the understanding of the spectral and 
dynamical properties of the associated random Schrodinger operators has been an ongoing challenge 
both for mathematicians (monographs include 043) an d for physicists (see for recent reviews). 
A nice geometry to study this question is that of an infinite regular tree. More specifically, given 
such a regular rooted tree graph Tk of fixed branching number K (each vertex has K+l neighbours, 
except the root which has K neighbours), the Anderson model on this graph corresponds to the 
operator acting on the Hilbert space £ 2 (Tk) 

H K>t {Lj) = -tT K + V(u) (1) 

where Tk is the adjacency matrix of Tk, t > is the hopping strength, and the "disorder" V is a 
random potential, i.e. a set of random variables V x indexed by the vertices of Tk- 

The issue of localization is related to the nature of the spectral measures ^K,t,5 x {dE, oj) associated 
with HK,t{u) and the Kronecker function 5 X € £ 2 {Tk) localized at x € Tk- More precisely, the 
different spectra of HK.t{^>) are associated with the components of the Lebesgue decomposition of 
the measures /j,K,t,s 3 ,{dE,uj): pure point, singular continuous or absolutely continuous. Ergodicity 
and equivalence of the local measures [7| ensures that the supports of the latter are almost surely 
non-random 043 and do not depend on x G Tk- The operator Hkj{u) exhibits spectral localization 
in an interval I C M if the spectral measures HK,t,S,c {dE, uj) associated with x € Tk are almost 
surely all of only pure point type in I. A stronger statement [3,11] is to say that HK,t{w) exhibits 
exponential dynamical localization in /: roughly speaking, this means that a state initially localized 
in space and in energy (within /) , and evolving according to the Schrodinger dynamics generated by 
HK.t{u), has a probability of being found a time r later at a distance R that decays exponentially 
fast with R, uniformly in r (we will not need this in the following, but a more formal definition 
can for instance be found in [9(). An opposite behaviour is absolutely continuous spectrum, which 
is equivalent to the statement that current injected through a wire at the corresponding energy on 
a site of the lattice will be conducted through the graph to infinity El ■ In the particular case 
of the regular tree, this transport happens ballistically [HI ]. 

The first analysis of Anderson localization on the infinite tree was performed by Abou-Chacra, 
Anderson and Thouless [l3l [l4| , later followed by many works from the physics community, see 
[ToL [H|-im and references therein. Besides the opposition between extended and localized spec- 
trum, recent topics of interest concern the nature of the extended and localized phases and of the 
derealization transition. From a mathematical point of view, the main results include analyticity 
of the density of states for random potentials close to the Cauchy distribution [Uj], spectral and 
dynamical localization of the spectrum for sufficiently strong disorder (2f| , and the persistence 
of extended states at sufficiently low disorder [IJ [28|, [29| . More recently, a major progress has been 
made by Aizenman and Warzel [9] by providing a criterion for the presence of extended states, 
which is (almost) complementary to the criterion for localization previousl y fo und in [2rl l27l ] . This 
allows them to work out the small disorder behaviour to a greater extent |30l Hl[ . 

An interesting limit in which the phase diagram of the Anderson model can be computed more 
explicitly is the large connectivity one K — > oo. This case has been considered early, see the 1974's 
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paper by Thouless [IH for a review of the different predictions available at that time. However, 
since then, no progress had been made and this large K limit was still considered unresolved 
[2l], [30l - The main contribution of this paper is to confirm the result obtained by Abou-Chacra, 
Anderson and Thouless for particular random potentials in 133, but on a rigorous basis and 
for generic disorder distributions. It is worth pointing out that this result can also be derived 
from the "best estimate" of Anderson's original paper Our proof is essentially based on the 
criterion for extended states given in [jj. Let us finally note that this large connectivity limit is 
not only of academic interest, since it has recently been argued to give a good approximation of 
the localization transition of many-particles systems (the connectivity being related to the number 
of elementary excitations, and the regular tree an approximation of the Fock space of the system) 

MM. 



We now turn to the presentation of the result derived in this article. We will assume that the 
random variables V x are independent and identically distributed, with a probability density p that 
satisfies: 

(A) p is piecewise Lipschitz continuous, meaning that there exist N p £ N*, —00 < eti < b\ < 0,2 < 

62 ■ • • < a/Vp < fr/Vp < 00 1 an( i a constant C\ > such that: 

- supp(p) c Uj(aj,6j) 

- Vze [l,Ag,V(«X)e { ai M) 2 Ap{v) - p{v')\ <Ci|«-«'| 

(B) There exist <; £ (0, 1) and C s > such that for all »eR*: 



C (2) 



Under Assumption (A), it follows that <^ could be taken as close to as desired; but we shall 
not use this in the following. 

(C) If supp(p) is unbounded, then for all k < 00, inf |„|<fc p{v) > 0. 

These conditions are for instance satisfied by linear combinations of Gaussian, Cauchy or piecewise 
constant distributions. It can easily be seen that they imply the conditions A to E of [9] (our 
last condition is stated only for this purpose and we will not use it in the following). Standard 
ergodicity results also imply that the spectrum of HK,t{u) is almost surely given by the set 
sum of the support of V and of that of —tTx '■ [—2t^/K, c lt\f~K]. 

For the sake of simplicity and to state an appealing result, we first state here our theorem under 
an additional assumption that we do not prove, but which is natural to expect, and was already 
suggested in Q. Since the statement of the latter requires the definition of additional quantities, 
we defer it to the beginning of next section, along with a more conservative statement of what we 
actually proved. 

Theorem 1. Assume that tjt = K £ K , that E £ R is such that p{E) > 0, and that Assumption 
1 holds true. Then there exists a constant g c (E) such that : 

— if g > g c {E), there exists Se > such that for K large enough the spectrum of HK,t K {^>) is 
almost surely absolutely continuous in the range [E — 5e,E + Se]; 

~ if g < gc{E), there exists 5e > such that for K large enough the spectrum of Hx,t K (^) is 
almost surely dynamically localized in the range [E — 5e, E + Se]- 



- ME) (3) 



Moreover g c {E) is given by 



Note that in this regime where the hopping term t^ is asymptotically much smaller than l/y/K, 
the density of states of the model converges in the K — > 00 limit towards the density of disorder 
p. Also note that, as observed in Q, this result also applies to the corresponding operator on the 
fully regular tree graph (or Bcthc lattice) in which every vertex has K + 1 neighbours. 
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The plan of our paper is as follows: in the next section we first state and dicuss in more details 
the result that we obtain before writing down an explicit expression for the "free energy" function 
introduced in @. In Section |m] we derive upper and lower bounds on this quantity that match 
when the connectivity K + 1 becomes large, allowing to prove our result. Section IIVI presents a 
numerical confirmation of our result, while we draw our conclusions in Sec.[V] Some more technical 
computations are deferred to an appendix. 



II. THE "FREE ENERGY" FUNCTION ip K ,t{s,E) AND ITS COMPUTATION IN THE 

LOCALIZED PHASE 



A. Statement of our result 

Following [i| [2|| [35[ , we will be interested in the asymptotic behaviour of the Green's function 
between the root and a site at distance L that we denote L; the latter is defined by: 

G K , t (0, L,E + i Vt w) = (S , {H K ,M -E- iv^Sl) (4) 
where r\ £ R!j_. Introducing, for s < 1 

VK , t{ s, E) = hm hm l° g ^ t (0,L, g + ^)|* 

rfSfi L 

where E denotes the average with respect to ui and the limits are well defined for almost all E, and 
commute and ipx,t{^, E) = lim S/ *i fK,t(s,E) (the limit is again well defined, see [jj), we can 
state a weak version of the main r 

esultsof|a|2fl|2Z|: 



Theorem 2 (SUSHI), a ssume that p satisfies assumptions (A-C). Let I C R be a fixed interval. 
Then: 

— if for Lebesgue almost all E 6 I, ipK,t(^,E) > — log if, the spectrum of Hx,t{w) is almost 
surely absolutely continuous in I, 

— if for Lebesgue almost all E G I, <fiK,t (1> E) < — log if, the spectrum of Hx,t(w) is almost 
surely of pure point type in I. Moreover if ess sup Se/ fK,tO-i E) < —logK, then i?x,t(w) 
almost surely exhibits dynamical localization in L . 

Our result is then as follows: 

Theorem 3. Assume that tx = xiogK an< ^ ^at E € R is such that p(E) > 0. Then there exists 
a constant g c {E), such that for all K at which the operator HK,t K {u) exhibits spectral localization 
(i.e. pure point spectrum) in a neighbourhood V of E : 

— */ 9 > 9c(E) and K is large enough, there exists 5e,k > such that tpK,tO-t E') > — log A" 
for almost all E' £ [E — Se,k, E + Se,k] 

— if g < g c {E) and K is large enough, there exists Se.k > such that 
ess8up E , e[E _ 5EKtE+5EJ(] tp K> t(l,E') < -log IT. 

Moreover g c {E) is given by (0). 

Intuitively, the first part of the theorem can be thought of as a proof by contradiction of der- 
ealization, while the second part shows that localization is "self-consistent". The reason for this 
complicated formulation is that we will compute tpK,t{s, E) conditioned on the fact that E lies in 
the localized phase. However, if we make the following additional assumption on ipK,t{^-t E), we 
can bypass this difficulty to obtain the nicer Theorem [T] 

Assumption 1. Assume that for all K > 2 the function (t,E) — > ipK,t(^,E) is continuous in 
(t, E), and that for all t, E <pjc,t{^-j E) + log K has only isolated zeros. 
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Let us explain how we derive Theorem [T] from our Theorem [3] and this additional assumption. 
We fix S E > such that g > sup E , e[ E-s E ,E+8 E ] g c (E'). For t K = K ^ g K and A large enough, 
one cannot have ipK,t{^-,E') < — logA" for E' € [E — S E ,E + S E ] since this would imply, using 
the continuity of E — > y>ir,t(l, A), Theorem [5] and then Theorem^ that tpK t t(l,E') > —logK. 
Hence the only possibility is that ifK t t(l,E') > — log A'. Using the second part of Assumption 1, 
it follows that for A large enough, (pic,tO-i •) > — log A" almost everywhere in some neighbourhood 
of E' . Using again Theorem [2j we get the first part of Theorem [TJ 

For the second part of the theorem, consider for S E > the set L e,s e — {t € R,Vt' < t,VE' € 
[E — 5 El E + 5 E ],ipK,t(^,E') < —log A"}. For t small enough (strong disorder regime), the 
spectrum is known to be of pure point type [HI, H3j- Thus it follows from Theorem [3] that L E _$ E 
is not empty. Using the uniform continuity of t — > <pK,t{l, E) on every compact subset, it follows 
from Theorem [5] and [3] that sup L Et g E is the smallest t such that, assuming that HK,t exhibits 
spectral localization in [E — S E , E + S E ], swp e ,^ e _ SeE+Se -i tpK,tO-,E) — — logA. Taking S E small 
enough and A large enough, this ends the proof. 

Finally, let us note that the bounds that we shall obtain on <pK,t(l,E) could in principle be 
used to deduce lower and upper bounds on a mobility edge (a curve in the (t, E) plane separating 
localized from extended states) also for finite A (under Assumption 1); however, without further 
work to optimize them, the latter are very far away from one another unless A becomes extremely 
large, and therefore we did not pursue this direction. 



B. Computation of ifK,t(s,E) in the localized phase 

To determine the boundary of the localized region, it is enough to compute the function (px,t(s, E) 
for s close enough to 1. To do this, let us first explain how one can compute the Green's function 
((31) between the sites and L. We denote 1, 2 . . . , L — 1 the sites on the (unique) path between 
and L, and L + l one of the neighbours of site L different from L — 1. Given a site x, the set of its 
neighbours y such that dist(0,y) = dist(0, x) + 1 will be denoted N£ . We will use the following 
expression for Gx,t(0, L,E + ir],u)) (see for example [§]): 



G K , t (p,L,E+iri,w) = 7 I] 



t 



t *J- Q V x (uj) - E-i-q-t 2 J2 y eM+ ^K,t(y,E + i<q,u) - t 2 T K . t (x + 1, E + iri,uj) 

(6) 

This expression is standard and can be derived using resolvent identities. In this equation 
TK,t(y,E + ir),uj) is the Green's function on site y and at energy E + ir\ for the operator H^.t 
obtained from Hkj by removing the path between the root and y. It satisfies the recursion 
equation d El HI HI] : 

r K , t (y, E + ir],w) = — — — ; -—^ — - — — — (7) 

Vy{u) — E — ir) — t 2 £ ze A# rx, f (z, E + iri,u) 

This is a random quantity because of the randomness in the V y . If we introduce the probability 
distribution PK,t,E+iri of TK,t(y, E + irj, lo) with respect to the underlying probability measure, we 
arrive at the following recursive distributional equation: 

r = k — (8) 

V - E - ir) — t 2 Y,?=i 

In this equation, both T on the left hand side and the Ti on the right hand side are independent 
random variables drawn from pK.t,E+ir), while V is distributed with density p. This equation is 
known to admit a unique solution as long as r\ > [36| . 

In the localized phase (i.e. if the spectrum is almost surely pure point in some neighbourhood 
V of E), for almost all realizations of the disorder, lim^o Tic,t(y, E' + ir], to) exists and is real for 
all sites y and almost all E' £ V. Conversely, for almost all E within the localized phase, and 
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almost all realizations of the disorder, lim^o ^K,t(y, E + 177,0;) exists and is real for all sites y. . 
Therefore, for an energy E such that the above property holds, the distribution of Tfc,t(y> E+irj, uj) 
becomes supported on the real axis, and can be replaced by: 

r = X — R — (9) 

where T is now real, and the solution of this equation has to be selected as the limit of the one of 
([8|) when r\ \ 0. Since the random variable V admits a density, its sum with any other random 
variable also does [38[ , and we deduce from © that T admits a probability density, that we denote 

PK,t,E- 

Because s < 1 the expectation value of \G(0,L,E + m],uj)\ s conditioned on the values of the 
disorder on all the sites x 4 jO, 1 , . . . , L} is uniformly bounded with respect to this disorder (this 
argument can be found in [9, 26]) 

E [\G Ktt (0, L,E + iri,^)\ s \{Vx}x^{o,...,L}] < (^j^) + t sL (10) 

This can be obtained here by noting that only the first x factors in the product defining © depend 
on V x . (|10p then follows from averaging the absolute value of © raised to the power s over the V x 
from x — up to x = L, and using repeatedly the bound E [l/\V — a\ s ] < 2 s (cf. Eq. (A5) of 
[§]). It follows that, using the convergence in distribution of TK,t(y, E + irf) towards TK,t{y, E), we 
can substitute TK,t{y,E) for TK,t(y,E + ir/) in \im T] \ i o'E\GK,t(0, L, E + i?7,u;)| s . Introducing the 

probability density PK,t,E of V — E — t 2 ^2i=i 1 (with the Ti independently distributed according 
to the solution pK,t,E of ©), we arrive at the expression that will be our starting point for the 
following study: 



limElGA^O^i^r^)! 8 = i / TT 



L 

I' 

PK,t,E(e x )de x 



(11) 



\e x -t 2 r Ktt (x + l,E)\s 

PK,t,E( r K,t{ L + 1, E))dT K>t (L + 1, E) 

where we recall that, with these notations, 

IV,t(z, E) = 1 (12) 

e x - t 2 T K . t (x + l,E) 

In particular TK,t(x,E) depends on (and only on) all the "energies" (e y ) x <y<L, and on Tx,t{L + 
1,E). We now introduce 



— > J X oc fcK,t,s,E(x,y)a(y)dy 



(13) 



with 



t 2 - s ( t 2 

ICK,t,s,E(x,y) = 1 .„_ OK,t,E V ) ■ (14) 



\x\ 2 s \ X 



Then it follows from (JXTJ) and changes of variables e x — > - _ t 2 Tli t ( x+ i m that: 

limE\G K A0, L,E + z Vi uj)\ s = h\F^ E (p K , t , E (- ' A 2 ))lli (15) 



1 Here it would be tempting to say that this property holds for all E in the localized phase, but we do not have a 
proof of this statement. One can however show that any E in the localized phase is almost surely not an eigenvalue 
of Hj(,t(u)', this follows from the simplicity of the spectrum within the pure point phase [37l | and ergodicity 
in particular Prop. V.2.8). 
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where the initial condition is given in terms of the probability density px,t,E of the solution F of 
@. Finally, we arrive at: 

^(l.-E) = lim lim log (\\Fk it , BiE (pK,tA- ' /* 2 ))llf~) (16) 

which, we recall, is valid for all E such that lim^o Tx,t{y, E + if], uj) exists and is real for all sites 
y and almost all uj. 



C. Connection with the criterion of Abou-Chacra, Anderson and Thouless and sketch of 

proof 

In the next section, we shall compute the sign of ip K « (1 , E) + log if when g is fixed, 
K is large, and under the assumption that Eq. p^|) holds, which will be our Fheorem [3J But 
before that, let us briefly explain heuristically how we can recover the critical condition given by 
Abou-Chacra, Anderson and Thouless in their paper [l3j from this point. 

Let us proceed as if the operator Fx,t,s,E was a finite dimensional one. Then using the positivity 
of the kernel K,K,t,s,E, we could apply the Perron- Frobenius theorem and conclude that ipK,t(s, E) = 
\og\K,t,s,E where \K,t,s,E is the largest eigenvalue of the operator Fx,t,s,E- The critical condition 
would then correspond to \ic,t,i,E = V-^' ana - ^ ne critical value of the hopping would be given by 
the largest value of t for which the equation 

Kt 2 - S f ( t 2 \ 

a ( x ) = u |2-s J P K ' t > E [~ y ~ ~) a (y) d y ( 17 ) 

possesses a solution up to s = 1. Making the substitutions t — > V, pkj,e{z) — > Q(—z) and 
s — > 2(1 — /3), we obtain Eq. (6.7) of 13]. A closer look at our derivation of ()16|) reveals that 
the eigenvector of largest eigenvalue of (|17p . a(x), has a simple interpretation: it is, up to a 
normalization, the limiting conditional expectation value of \G(0,L,E + i?7,u;)| s , conditioned on 
TkjO-, E) being equal to —x/t 2 . A very similar interpretation can be given for the dominant 
eigenvector of the adjoint equation of fp~T|) . used in [13]. Note that the interpretation that we 
obtain for these eigenvectors is somehow more explicit than the one given in [l3j . where they 
appear as being related to the amplitude of the tails of the distribution of the local Green's 
function G(0, 0, E + if], uj), assuming a power-law tail with an exponant 1 + ft- 

Here it is worth showing the numerically determined eigenvector of largest eigenvalue of a finite 
dimensional approximation of K,x,t,s,E, close to the transition (see Fig. [1]). The latter is very well 
fitted by a(x) oc l/\x\ in a wide range of x, except for some cutoff near x = and x — > ±oo. The 
general idea of our proof is then simple: if the operator Fx,t,s,E was acting on a real Hilbert space 
of finite dimension, then we could use the Collatz-Wielandt formula [39] to derive 

. , / JC K ,t,s,E(x,y)v(y)dy . J K. K ,t,s,E{x,y)u(y)dy 
sup ml -—^ = \K,t,s.E = ml sup -— (18) 

vhO v(x X )>0 ^° x U{x) 

Here the order is the partial order defined by u -< v iff u(x) < v(x) for (almost) all x. Then, an 
asymptotic study of the conditions (fT8| lor K large, using well-chosen test vectors, would allow us to 
conclude. However, Fpc.t,s,E acts on a Banach space ol infinite dimensiomand is not compact (il it 
was, we could use natural extensions ol finite dimension spectral results [401] ). We shall circumvent 
this difficulty in the next section by using an elementary version of the Collatz-Wielandt formula. 
Note that the fact that px,t,E appears in K,K,t,s,ei instead of the bare density p, should not be 
relevant in the large K limit, as was already noted in the early works [III, UH - even though it will 
add some technical difficulties in the following. 

Finally, our method is not (in its present form) efficient for s < 1, but one can still study Eq. 
(fTT|) numerically to estimate the free energy function tpK,t(s, E). What one finds is that ipx,t{s, E) 
is equivalent to s log t for s going to (as was noticed in [31j for the Cauchy disorder case) . Such 
a linear approximation to ipK,t(s, E) would lead to a critical value of the hopping scaling like 1/ K. 
However, lor s > 0, there are subleading (and disorder dependant) corrections to this linear form 
that go as log(— \ogt) for small t, giving rise to the correct scaling tj< — g c /(K log K) for the 
critical value of the hopping strength. 
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FIG. 1: jx|a(x), where a(x) is the eigenvector of largest eigenvalue for a discretization of the kernel K,K,t,s,E, 
for K = 2, s = 1, E = 0, and two densities of disorder p: a uniform distribution p(v) — jl(t 6 [—1, 1]) 
shown with a solid black line (for t = 0.11), and a standard Cauchy distribution p(v) — w ^ v i +1 ^ shown 
with a dotted red line (for t = 0.23). Right: detail around x = 0. a(x) has a finite limit when x — > 0. In 
each case t is chosen close to its critical value (the precise value of t has no influence on the global shape of 
the eigenvector). a(x) was obtained by iterations of the kernel K,K,t,s,E, using a discretization of the real 
line in 16000 elementary intervals. 



III. ASYMPTOTICS FOR ip K>t (s,E) IN THE LOCALIZED PHASE 

As explained above, our proof is essentially based on the following elementary proposition: 

Proposition 1. Let (X, || • ||) be a normed vector space with a partial order ^ compatible with the 
product by a non-negative real (u^v and A > =>■ Ait r< Xv), and such that 

Wu,veX, (O^u^v) ||u|| < \\v\\ (19) 

Let F be a linear application from X to itself preserving ^ (u -< v =>■ F(u) -< F(v) ). Let u and v 
be non-negative vectors (0 ^ it,0 -< v), different from the null vector, satisfying respectively: 

F(u) -< A u 

b (v) >z (JL V 

for some A, fi > 0. Then for all vector p£X such that there exist a, b > for which 

av <p<bu , (21) 

it holds that: 

ix < liminf ||F n (p)|| 1/n < limsup \\F n {p)\\ 1/n < A (22) 

n ^oo n — 

The proposition easily follows from ([20]) and (|21|) by successive applications of the operator F. 

In the following, without loss of generality we prove Theorem [3] for E = - the generic case can 
always be recovered by considering the shifted density p(- + E), which satisfies the same hypotheses 
(A-C) as p. We will assume that E — lies in the localized phase, meaning that (K, t) <G Loc 
where 

Loc = {(K, t), there exists a neighbourhood V of E = such that the spectrum of FLk.x 

exhibits spectral localization in V.} 

Because of the technical reasons mentioned in Sec. Ill Bt we still have to consider a small interval 
of energies around E = 0. To do so, let us first consider io such that € (ai Q ,bi ) (recall that the 
{fli, b{) were defined by assumption (A) on p as intervals of R on which p is Lipschitz continuous) 
and introduce m £ (0, 1/2) such that [—to, to] C (aj Q ,&j ). For (K,t) £ Loc we also fix < 6j(,t < 



s 



min(l/K, to/2) such that the spectrum of Hpc.t is localized within [— Skj, #K,t]) and we finally 
define 

Ek t — {E G [—5k t,5K t], lirn t {y, E + it], uj) exists and is real for all site y and almost all uj} 

v\o 

(24) 

This allows us to use formula (|16j) to compute ipK.t{s,E) for E € fx,*, while the bounds on 
and to are unessential and needed only for technical reasons. As already mentioned, we know that 
[—5K,ti 5kj] \ &K,t has zero Lebesgue measure (and we conjecture that this is actually an empty 
set). 

Finally, we also fix G > such that, for all K, the hopping t belongs to the interval [0, G/K]. 




FIG. 2: Sketch of the various quantities defined in the text for fixed (K,t) £ Loc: a,i a ,bi are defined by 
condition (A) on p and £ {a io , b io ); m £ (0, 1/2) is such that [— m, m] C (a io , b io ). The range of energies 
I such that the spectrum of Hkj is localized within / is shown by a red line, and 5k, t £ (0, 1/K) is then 
chosen such that [— Sn.t, &K,t] C [— m/2,m/2]. Finally £K,t is the subset of full measure of [—Sx.t, 5K,t] 
defined by |gg) . 



We will now apply the proposition to the operator FK,t,s,E on the normed vector space L^R), 
with the partial order ^ defined by u ■< v if and only if u{x) < u(cc) for almost all a; G M. Since 
FK,t,s,E is linear and has a non-negative kernel, it is an easy check that it indeed preserves the 
order ■<. 



A. Lower bound on ifK,t{s,E) 

We consider for (K,t) G Loc and E G £-K,t fixed, the following "test-vector" v, for a given 
A G (0,to/4): 

= if f<N<i (25) 

w \ otherwise v ' 

As discussed in Sec. Ill C| for A well-chosen, this vector is very close to a numerically determined 
dominant eigenvector for |a;| G (0, 1); it will turn out that the expected shape of the eigenvector 
for | a; | > 1 is not needed for this variational computation. 

We shall prove in the appendix (cf Eq. (|A8[1 ) that the vector Pif,t,s(— T/t 2 ) is uniformly bounded 
away from zero on any interval of E*. Since v(x) has bounded support and is bounded, (l2"Tj) will 
then be satisfied for some a > 0. Moreover, v(x) satisfies FK,t,s,E{v) h l^v as soon as: 

. nf J lCK t , s M X[ yHy)dy 

|x|e[A,i] v(x) 
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Hence 



inf 

|x|e[A,i 



IK 



K,t,s,E 



(x,y)v{y)dy 



We compute, for x G [A, 1]: 

J K, K ,t,s,E(x,y)v{y)dy 
v(x) 



v{x) 



PK,t,E 



A<\y\<l 



-y — -n d y 
\y\ 



(27) 



(28) 



The idea is now to take A small when K — > oo, so that the above integral becomes dominated 
by what happens for y small. Since moreover px,t,E should resemble p(- + E) for K large and t 
accordingly small, the integral in (|28|) should become equivalent to —2p(E — t 2 /x)\ogA. It will 
be more convenient to proceed in two steps, and to first substitute px,t,E for px,t,E, where pK,t,E 
is the probability density of V — E — t 2 YliLi 1\ (there is K terms in the sum instead of (K — 1) 
for pK,t,E, and the Ti are drawn independently with density PK,t,E)- This is made possible by the 
following lemma, whose proof is deferred to the appendix: 

Lemma 1. There exist two constants C\,Ci > such that for all K > 2, t < Gj K ,E G £n.t and 
z G K: 

-2 (p K ,tA*) - -7= ] log A-2C 2 < / PK .t,E {-y + z) ±-dy < -2 (p K ,t,E(z) + 4= J log A+2C 2 

V VKj JA<\ y \<i \y\ V vkJ 

(29) 



Replacing in (|28|). we obtain that: 
fK 



hJ..-,Eix,y)v{y)dy Jf 



v(x) 



-2 PK, 



,t,E 



This leads, using (j2"T)) . taking s /*• 1 and setting A = i 2 /a, with a G (0, 1), to: 



tp K ,t{l,E) > log t + log 



~4 m K:t (a) - 



\[K 



logi + 2 M K , t {a) 



Ci 
K 



log a - 2C 2 



where = max(:r, 0) and we introduced 



mf p K ,t,E(z) 

\z\<a 
E£S K ,t 



M K>t (a) = sup pk.i,e{z) 

\z\<a 

Ee£ K ,t 



(30) 



(31) 



(32) 



Remembering that pK,t,E is the density of 1/ — E — t 2 ^2%=i ^ii i n the limit of K large but with 
i < G/K it should get close to the bare density p(- + E). More formally, we prove in the appendix 
(cf Eq. (|A25[) ) that, in the limit of K large, with tx depending on K and satisfying tx < G/K, 
and a < m/2 2 : 

(33) 
(34) 
when 



(35) 



lim niK.t K (ex) = m(a) = inf p(z) 

if— >-oo " M< a 

lim Mr-^ (a) = M(a) = sup p(z) 



K->oo 



\z I <a 



From (|3"Tj) and Q33H34p we deduce the asymptotic behaviour of ipK,t K (l, E) near E 
Ik = and g is large enough: 



,<Pk,- 



.(l,E) + ]ogK > 



if log is" 

Vg > -. — -, lim inf ess inf 

4m(a) if^oo Ee[-5jc,t K ,5K,t K j " " 

(^KTCiTr)e Loc 

Since this is valid for a arbitrarily small, we deduce the second part of Theorem [3J 

Vo > - — 7-7, lim inf ess inf w K 

4p(0)' *-►<» BeM 6 J*' A >*i 

(^irnfiK)6 L ° c 



,(1,J5) + log A' >0 



(36) 



One may have to take K sufficiently large or G sufficiently small to satisfy the constraints A < m/4 and a < m 
at the same time; for the sake of clarity we keep this condition understood in the following. 
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B. Upper bound on ifK,t(s,E) 



To derive an upper bound on <fK.t{s, E), we consider for A £ (0, m/4) and A > the following 
vector u e L 1 (M): 



u{x) 



2-s 



x 

A 1 



if |a;| > 1 

if A < |af| < 1 

if Id < A 



(37) 



(<r was defined by condition (B) on p). Using the bounds on pK,t,E derived in appendix (Eq. 
(|A2IIA7P ). u satisfies ([21]) with respect to PK,t,E{— • /t 2 ) for some b > 0. Then we have, for \x\ < A 
(and denoting as before A = t 2 /a): 



J JC K ,t,s,E(x,y)u(y)dy _ t 2 



u(x) 



ll-s+S 







/ PK,t,E ( 


- y -x)- 


yi»i<A v 





\v\ 



1 j A? 



A<|y|<l 



l 

¥\ 



A ? 



i»i>i 



< f 



A s 



£2-2s+2? 



sup 



l-a+S 



2 ( - - ilogA ] sup (y + z) |z 

? ^ / |»|<i 



ll-s+S 



PK,t,E {y + z)-r-. 

y \>i \y\ 



2-s 



-dy 



(38) 



We now use the following lemma, proved in the appendix: 

Lemma 2. There existC 3 ,C 4 € R swcft ifcat, /or a// K > 2,t < G/K, E e an-d s > (l + ?)/2; 

(39) 



sup sup pjif.t,^ (y + z) jz 1 s+? | < C 3 
2 |y|<i 

sup / PK.t.B (y + z) dy < C 4 

* j\v\>i \y\ 



This gives by replacing in (|38j) . and for s > (1 + ?)/2: 



J K. K ,t, s ,E{x,y)u(y)dy 
sup t— r <ia 

|s|<A «W 



2 4, 2, \ „ 1 

-- -logt+-loga) C 3 + -C 4 



(40) 



(41) 



On the other hand, for \x\ > A, and using again Lemma [T] 

±2-s 



J K. K ,t, s ,E(x,y)u(y)dy t 2 



u{x) 



|A| 



l-s 



2— \\pK,t,E | 



A<\y\<l 
t 



2A M^_ 2 (p K , t , E [ 

? V V x 



1 2/1 

logA + 2C 2 + l 



(42) 



To go from the first line to the second, we used that ||pK",t,£;||oo < IIpIIoo! this follows from the 
fact that PK,t,E is the density of a sum of random variables that contains a random variable 
V distributed with density p, and from the bound on the convolution of two densities /i,/ 2 : 
ll/i * /alloc < II /l ||oo- Taking A = 2C 3 ^, and keeping t<G/K,$g$ will be larger than gT) for 
K large enough (independently of s > (1 + ?)/2). Therefore, for if large enough: 



<p Kr t(l,E) < logt+log 



^ "(0) IIP | lo ° " 4 ( Mif - t(a) + 7f ) log * + 2 (^*( a ) + 4f ) lo s a + 2C 2 

(43) 
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Using again (|33H34|) . we deduce the asymptotic upper bound (for a < m): 



4M («) , K~>oo Be[ - 5 S ] »'■-» 



Since this is valid for a arbitrarily small, we deduce the first part of Theorem [3] 

V # < 7~77^' limsup esssup (p K .^_^(l, E) + log K < (45) 

M u ) , Ee[-<5 K , t „,« K , t „] 



IV. NUMERICAL RESULTS 



In order to demonstrate the agreement of our results with numerical studies, we performed 
several numerical simulations. To simplify the discussion, we will assume in this section that there 
exists, for the energy E considered, and for every connectivity K + 1, a unique number g c (E,K) 
such that the spectrum of H K _s_ is localized in a neighbourhood of E if g < g c (E,K) and 
absolutely continuous in a neighbourhood of E otherwise. In particular g c {E) = limx^-oo 9c(E, K). 
We will focus on results obtained for two choices of the probability density p of the disorder: 
the uniform distribution p(v) = ^l(i> € [—1,1]) and the Cauchy distribution p(v) = ^njk^ ■ In 
fact, in both cases, as we will explain later, we have a prediction not only for the leading order 
of (fiK,t(l,E) within the localized phase, but also for the first subleading correction. Moreover, 
in the Cauchy case, it is well-known that the density px,t,E can be computed exactly, which will 
simplify some computations. 

The simplest quantity to determine numerically is the largest eigenvalue of a finite dimensional 
approximation of the kernel K-K,t.s,E, with the additional assumption that one may substitute p for 
PK,t,E in K-K,t,s,E- Since this (modified) kernel is positive and irreducible, the latter can be found 
by iterations - a few are enough. One then has to let the number of discretization steps go to 
infinity in order to "correctly" approximate the original continuous operator (this remains however 
uncontrolled). Since we do not use the exact density of states px,t,E in this case, this is still an 
approximation of the operator Fx,t,s,E which is expected to give good results for the mobility edge 
g c {K, E) only as K becomes large. More precisely, one expects the finite K corrections due to the 
use of p instead of pK,t,E to be quite small (of order Kt 2 — 0{l/{K\og 2 K)) [H,[IE|]); in particular 
much smaller than those in powers of 1/ log K that we expect to find regardless of the precise form 

of p. We denote gi A \K,E) the critical value computed with this procedure. 

In order to deal with the true kernel K,K,t,s,E, one has to determine the modified density of 
states pK,t,E first. For the Cauchy distribution this can be done analytically [id. 1X31] : in general 
one first has to solve the recursive equation ^ with some population dynamics approximation 
[ToL Il3l . Il9l . [2lj - the precision of the latter resolution is not critical in this case, and moderate 
population sizes and number of iterations are enough to get a good approximation of PK,t,E- 
Plugging this into K,K,t,s,E and following the same route as before to find the largest eigenvalue 

of this kernel this gives a critical value gc (K, E) which is exact (assuming infinite numerical 
precision) and is the most convenient to compute. It can be seen in Table HI HI] that this value is, 

as anticipated in the previous paragraph, very well approximated by (K, E) even for moderate 
values of K. 

A different way to proceed is to look for the localization transition directly on the tree. In this 
case it is actually much more convenient numerically to consider the "quenched free energy" 



4> K ,t(s,E) = lim lim —Elog 

rf\0 R~>oo R 

(note that compared to @,the log is outside of the sum). The latter can be computed, for any 
s, with the cavity method [23]. The localization transition was usually found [13 EH O, H3 by 



E 

xeT K 

dist(0,x)=_R 



G K , t {Q,x,E + ir),u)\ s 



\ogK 



(46) 
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looking at the critical value of t such that <j)K,t{2,E) equals — log if, because this corresponds 
to the apparition of a non- vanishing imaginary part in T solution of (151) when f) \ |, l23l . l3lj |. 
In this case, strong finite-size effects were mentioned [To|, EH [H, HH, HH and the values obtained 
were hard to reconcile with those obtained with the method presented above [HI]. However, at 
the localization transition, (j>K,t(s, E) does not depend on s for s € [1, 2] and is equal to <pK,t{s, E) 
[23d . One can therefore compute 4>K.t{s,E) as a function of t for any s g [1,2] and check when 
it reaches the critical value — log if. It turns out that doing this simulation for s = 1 strongly 
reduces the finite-size effects. Moreover it also allows to take r/ = from the beginning (within 
the localized phase), reducing the number of parameters on which the numerical results depend. 
In addition to that, we performed two finite-size scaling analysis. To explain them, let us call 
4>k'^ (s, E) the approximated value of 4>K,t{s,E) computed with R iterations of a pool of size 
N (the procedure is analog to that of [23[). We first take the R — > 00 limit at fixed N, using a 

finite-size correction as 0^ (s, E) = ^Kti 8 ' E) + a.K,t(N)/R and R in the range [1000, 10000]. 

Then we take the population size N to infinity, assuming logarithmic corrections of the form 

4>^ N \s, E) = tfr^°°\s,E) + b K j/logN + CK jt /(log N) 2 . We used this form because it fits 

relatively well the corrections, and also because it was predicted in [l9j] for this particular model; 

of course such slow corrections impoverish the quality of the numerical fit. For the purpose of this 

work, we used population sizes between 10 4 and 4.10 6 , and also took an extra average over a few 

independent realizations of the R cavity iterations, in order to reduce fluctuations effects. All of 

( o 

this leads to a computationally heavy procedure. However it gives critical values g c (A, A) that 
are quite different from the ones previously computed with this method [HI, HH , but very close 

to the exact value g c (A, E) computed with the procedure above. In particular this resolves the 
discrepancy that was observed in [T3[ . 

Finally we want to compare these results with analytic asymptotic predictions for g c (K, E). The 

first one is a straightforward consequence of the bounds of Sec. IIII1 and defines g$P [K, E) as the 
smallest root of 

log 5 - log log A' + log [-ip{E) (log g- log(Alog A))] (47) 

This expansion can in principle be improved by computing the next subleading term in ipKj(l, E) 
for A large. In the particular case of E = 0, one should recover the prediction of [l3| (Eq. (7.6)), 
obtained through an approximate solution of the eigenvalue equation (|17p for the operator Fjc,t,s,E- 

This gives, both in the uniform and in the Cauchy case, a critical value g c (A, 0) defined as the 
smallest root of: 



log g — log log A + log 



^ (0) I 1 + 24 (log g -log ( AlogA)r 1 (1 ° g5 " W K **V) 



(48) 



Note that there is no reason to expect this expansion, that may seem universal at first sight (the 
same numerical factor 7r 2 /24 appearing in both cases), to remain valid for other disorder densities, 
or for £ / 0. Also note that in the uniform case, this last result can be recovered through a 
systematic expansion of lpka{s,E) in powers of 1/logA, starting from Eq. (fTTT) 41]. In any case, 
we expect that 3 : 

5 f) (A, 0) = g c ( A, 0) + O ( (7^)2 ) ^ E) °) = 9c(K, 0) + O (^53 ) (49) 

where the 0{. . . ) terms should be rather small, even at finite A. On the other hand we expect for 
the critical value at finite A, g c (K,E), a very slow convergence towards its limiting value: 

, (K , E) ^ (E) __LW + (_i_) (50) 



3 In this particular case only; for general densities p and energies E, we expect the first correction to to go as 
1/ log if. 
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Hence it is worth considering the subleading corrections that led to (K, E) . The agreement 

between these values (in particular <?c (K, 0)) and the best numerical value gc (K, 0) can be seen 
to be very good in the uniform case, and also satisfying in the Cauchy case, even though the finite 
K corrections are stronger. 



K 


<^ A) (a,o) 




^ c) (A,o) 




gi E) (K,o) 


2 


0.150 


0.153 


0.154 


0.154 


0.149 


3 


0.187 


0.188 


0.189 


0.194 


0.187 


4 


0.207 


0.208 


0.204 


0.213 


0.207 


5 


0.220 


0.220 


0.219 


0.225 


0.220 


6 


0.230 


0.231 


0.227 


0.234 


0.230 


8 


0.243 


0.243 




0.247 


0.243 


12 


0.261 


0.261 




0.263 


0.260 



TABLE I: Critical values of the disorder for the uniform case and E = 0, for several numerical procedures 
(superscript A to C) and two asymptotic formula (superscript D and E), defined in the text. The numerical 
error for gi^ and g c B ^ is smaller than 0.001, while it is larger for g c C ^ (of the order of 0.005). In this case 
the asymptotic value is g c (E) = 1/2, and is approached extremely slowly. The connection with the model 
where V is uniformly distributed in [-W/2, W/2] and t = 1 is given by W c = — A" log if. 



A 


^ A) (A,0) 


ffi B) (A,0) 


^ C) (A,0) 




g?\K,Q) 


2 


0.317 


0.334 


0.334 




0.367 


3 


0.364 


0.372 


0.370 


0.418 


0.384 


4 


0.389 


0.394 


0.394 


0.423 


0.403 


5 


0.406 


0.410 


0.404 


0.432 


0.417 


6 


0.419 


0.421 


0.422 


0.440 


0.428 


8 


0.436 


0.437 




0.453 


0.444 


12 


0.456 


0.457 




0.470 


0.463 



TABLE II: Critical values of the disorder for the Cauchy case and E = 0, for several numerical procedures 
(superscript A to C) and two asymptotic formula (superscript D and E), defined in the text. The numerical 
error for g c A ^ and gc is smaller than 0.001, while it is of order 0.01 for gi ■ In this case the asymptotic 
value is g c (E) = 7r/4 ~ 0.785. The connection with the model where V has Cauchy distribution of 
parameter 7 and t = 1 is given by j c = — A" log A. 



V. CONCLUSION 

In this article we have computed the free energy function tpK,t introduced in H on tree graphs, 
in the large connectivity limit and within the pure point phase. Using additional unproved reg- 
ularities argument on ifK,t-, this gives the asymptotic scaling of the mobility edge in this regime. 
Interestingly, we have found that this rigorous approach gives back the criterion obtained in [l3[ 
using a self consistent equation on the tail of the distribution of the local Green's function. 

Our technique applies at any energy E at which the original density of disorder p(E) is continuous 
and strictly positive (and under additional assumptions on p that are probably mostly technical). 
An interesting extension of this work would be to understand what happens when p{E) = 0. This 
can happen in two ways: a possibility is that the energy E considered is such that p{E) — but 
some (one-sided) derivative of p does not vanish. Another appealing regime is when E is outside 
the support of the disorder density p. In this case, one may expect, following the bound given in 
f30j ] . a mobility edge scaling like 1/K in a range of energies at distance at most \j\f~K from the 
support of p. The behaviour of the free energy function tp K t and of the related localization length 
in this Lifshitz tail regime would be particularly interesting to understand. 

Besides the value of the mobility edge, another direction of investigation concerns the nature of 
the derealization transition on the Bcthc lattice. The difference between this transition and an 
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ergodicity breaking transition, related to more subtle properties of the spectrum such as eigenvalues 
spacings, was recently emphasized in [23| . It would be nice to obtain also an asymptotic for this 
transition when the connectivity becomes large; if this is to be done with the same approach, one 
would need to handle the complex solution to the recursive distributional equation. Finally, the 
connection with the finite dimensional Anderson model in the limit of large dimension (discussed 
for a related model in 42]) or on the hypercube are points that certainly deserve further studies. 
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Appendix A: Study of the real RDE 

This appendix is devoted to the derivation of various estimates for the probability density PK,t,E 
in the regime of K large and t accordingly small (t < G/K). Our starting point is the real recursive 
distributional equation ([9]), that we recall here for convenience: 

r = 1 ~ R — (Ai) 

In the following we shall always assume that this equation admits at least one solution, and speak 
about "the solution" of (|A1I) for any given solution of it. 4 Finally, recall that we always assume 
E G [-5 K ,t,8 K ,t\ vnthS Ktt < mm(l/K,m/2), m G (0,1/2) and [-m,m] C (a io ,b io ) 3 (cf Fig.E). 

We divide the results of this appendix in four (dependent) parts: first the study of the solution 
of equation (|A1[) . then its consequences on the regularity of px,t,E and the proof of lemma 1 and 
2, and finally a study of the convergence of pK,t,E towards p when K is large (and t small). 

We recall that we denoted PK,t,E the probability density of T solution of equation (|A1[) . Let 
us also denote p^t e ( res P- P^kTe) tne probability density of t 2 Y^f=i (resp. t 2 Ylf=i^i)- 
Finally, recall that we introduced the density pn,t,E (resp. pK,t,E 

) of V - E - t 2 Y*=i l (resp. 

V-E-t^f^T,). 

a. Probability distribution of the solution of the real recursive distributional equation 
In this first part we derive estimates onp^t.B that were needed to apply Proposition 1 in Sec. HID 



Since V admits a density, any sum of V with another random variable also does |3a | . In particular 
this holds for the denominator in (|Alj) . and thus also for the left hand side of (|Alj) . Hence pkj.e 
is indeed a density. Moreover it satisfies: 

< ^ (A2) 



PK,t,E(T) = J P ( kI e (x)p(~+E + x^j dx 



Using that, given n densities /i, /2, . . . , f n their convolution satisfies 

(fi*f2*--'*fn)(x)< sup /i(y)H h sup f n {y) (A3) 

y>x/n y>x/n 



4 This is not an issue since in the core of the text we use this equation only for energies E G £x t (defined in 
Sec. which guarantees existence of at least one solution to ||A1| |. while the choice of the relevant solution is 
made unambiguous by the unicity of the solution to the complex equation JS), as explained in Sec. Ill Bl 
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(and the same inequality with y > x/n replaced by y < x/n) we deduce that 

P g E (x)<^ 2 Mf (A4) 

From assumption (B) on p: Vt> ^ 0, p(v) < C^/\v\ 1+ ^ , we deduce the following bound for x > 1/2: 

C 2 1+? C 

sup p(x + E)< sup i, < — t—^ (A5) 

Ee[-m/2,m/3] Ee[-l/4,l/4] ( x + X + s 

Using this with (|A3|) gives that, for x > 1 

Pjr.t.E (z) < 4 (2 1+ ^C, + X 3 t 2 ||p||oo) 
Going back to V we obtain that for |T| < 1: 



(A6) 



PK,t,E(T) = ^f PK ,t.E (^j < 4 (2 1+ <C ? + X 3 t 2 ||p||oo) ^7 (A7) 

We also need the positivity of PK,t,E(T) to derive the lower bound of Sec. lIII Al To do this, first note 
that r — > PK,t,E(T) is a continuous function (except possibly in zero); this follows from the Lipschitz 
continuity of f>K,t.E, that will be proved below (cf Eq. (|A14[) ). Hence there exists [a_,a+] C R 
such that the restriction of PK,t,E to [a_,a+] is positive. Assume first that € [a_,a+]. Since 
there also exists 5 > such that the restriction of p{- + E) to [—6, S] is positive, it follows from one 
iteration of (|Alj) that V has a positive density near ±oo. Hence ^ I\j has a positive density on M 
(every real number can be written as the sum of K — 1 arbitrarily large real numbers), and with 
another application of (| Al|) . we deduce that the same holds for T. Using the uniform continuity of 
Px,t,E on every compact of R*, we deduce that 

Va,&eM,afr>0^ inf p K tE(^)>0 (A8) 
re[o,6] 

It remains to prove that one can assume G [a_,o: + ]. To do that, consider the successive images 
of a > by the map / : x -t e _ t l Kx - For e £ [-6,6], if a G supp(p K , /(a) € supp(p K .t, E )- 
Moreover, for e small enough, there exists n such that f 2n (a) € [— S/(t 2 K), 5/(t 2 K)]. Hence if 
a_ > 0, it follows from repeted applications of (|A1|) that T{ has a positive density in some 
subset of [—<5,(5]. The same reasoning can be done if a+ < 0, taking e < close enough to zero. 
Hence Y has a positive density near ±oo, and applying once again the recursive equation, we deduce 
that pK.t,E has a positive density near 0, which ends the proof. 



b. Proof of lemma 2 



We now turn to the proof of lemma 2 and 1. Recall that from now on it is assumed that 
t < G/K with G > fixed. 

It will be convenient to use the following union bound, for a < 0: 

X U • • • U | Tk^ < 



t 2 r * < a < 



ri < 



t 2 (K - 1) 



t 2 (K- 1) 



< (k - i)p r < 



t 2 (K-l) 



(A9) 



< 



G 2 ||p||c 



where we used, from (|A2|) . and for b < 



t(r< ^< i2 II^.G 2 ||p|| t 



< 



\b\ ~ K 2 \b\ 



(A10) 



1G 



From this, and with the help of (|A5j) . we deduce that for z > 1 
PK,t,E(z) = [ p(z + z' + E)p ( K-£ (z')dz' < 



[ p(z + z' + E) P <£-£{z')dz' < sup p(z + z' + E) + \\p\\ 00 P ( t 2 V Y % < -z/2 j 

J E6[-m/2,m/2] \ i=1 / 



< 



;'>-z/2 

, 2G 2 ||p|| 2 oc 



jl-K 2 

(All) 



Hence 



sup sup pxABfo + aOlsl 1- ** ^supp^B^J^+l) 1 -^ < 8 1+? G ? + 4G 2 || /9 || 2 <3 (A12) 

z>2|y|<l Z>1 

The same bound holds for z < —2, while for z £ [—2, 2] it is clear that 
sup UK 2 sup, Kl p K .t,E (y + z) z 1 ^ 3 ^ < ||p|| oc 2 1+<; (note that ||pK,t,B||oo < IIpIIoo! as explained 
in Sec. IIII Bl this follows from the fact that pK,t,E is the density of a sum of random variables 
that contains a random variable V distributed with density p, and from the bound on the convo- 
lution of two densities fi,f2- ||/i*/2||oo < ||/i||oo)- Hence we have proved the first part of lemma 2. 

For the second part we proceed in the same way: considering hrst \z\ > 2, we write: 

l-s-K 



PK,t,B(V + z) ^ 2 _ s dV = j ]y] > 1 PK,tAV + *) Z ^2-s d V + j M >! PK,tAV + Z ) Z \ y] 2-s d V 

\y\>z/2 " |y|<z/2 



1 1 y v k/2i i+? V2 y v 

y<M/2 

<2 2 -|zr 1 + 8 (2 1 ^Q- + G 2 ||HlL) prolog (1^1/2) 



(A13) 

where we used in the last line that lemma 2 assumes s > (1 + ?)/2. Hence z — » X| y | > i PK,t,E(y + 

z) ^yp-a dy is bounded uniformly in s, -ftT and t for |z| > 2. Since it is also clearly the case for 
\z\ < 2, we have proved the second part of lemma 2. 



c. Proof of lemma 1 

In order to derive lemma 1, hrst note that ||Px7e||oo < IIpIIoo; this follows from the same 
argument as the one given under (|A1~2|) . We now look at local regularity properties of pK,t,E and 
compute, for e— e' > and introducing B e y(p) — L) J ^ 1 (ai — E—e,ai — E—e')Ll(bi — E—e,bi — E—e'): 

\pK,t,E(e) - PK.t,E(e')\ < / \p{e + x + E) - p(e' + x + E)\p ( x} E (x)dx 



< 



W + x + E)- p{e' + x + E)]pf ) tE {x)dx+\\p\\ 00 / p ( p tl 



<Cx|e-e'| / P^izWdx + WpWl / dx 



< 



(C 1 + 2^,||p|| 2 00 )| £ 



Ci 

(A14) 

Hence pK,t,E is (globally) Lipschitz continuous, with a Lipschitz constant C\ that does not depend 
on K, t or E. 
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In order to deal with pK,t,E, we recall that E E [— m/2,m/2] where m was chosen such that 
[— m, m] C (a,i a ,bi ). Moreover, using the bound (|A4j) and t < G/K, we have that P^KtE( z ) — 
G 2 K\\p\\ 00 / z 2 . A computation similar to (|A14[) then leads to: 



Ve, e' G [-m/2, to/2] 2 , |/9jc a k(c) - PK,t,E^')\ < AK \ e ~ e 'l 



(A15) 



where Ael does not depend on K, t or E. Even though we expect that the factor K in the right 
hand side of (|A15j) could be dropped, the weak statement above will be enough for our purpose. 

We are now ready to prove Lemma 1. We first bound PK,t,E(e) with pK,t,E(e) as follows, for 
e e [-m/ A, m/4], and with the help of Eq. (|A10j) and (jA15|t : 



pK,t,E(e) = / PK,t,E(e + x)p K ,t,E(t 2 x)^- > 



dx 



> 



inf 



inf 



PK,t,e(e'W * 2 |r| < 



KVK 



PK,t,e(e') 1-2 



,G 2 |M| 



m\/ K 



Am o G 2 || /0 || 2 o 



(A16) 



i>K,.L(< )- sup _ pis.i.l(< ') + [j/»jj vt : [ / 2 |r ■ T7 ^=) < PM ( e) + 4^ + 2 G2|IHI - 



Hence, for A < m/4: 

f i \ dy f ~ i \dy 

/ PK,t,E[-y + Z)— PK,t,E(-y + z)—r 

Ja<\v\<i \y\ Ja<\ v \<i \y\ 



K mVK 
(A17) 



^UTO+^»)^(-logA) 



2||p|| 00 log(TO/4) 



(A18) 



Using that PK,t,E is Ci-Lipschitz continuous, we get: 



/ 

JA<\y\<l 

Lemma 1 follows from (|A18[) and (|A19|I . 



pK,t,E(~y + z)^-dy + 2p Kt , E {z) log A 



< 2Ci 



(A19) 



rf. Convergence of pK,t,E towards p for K large 



Finally, we shall be interested in estimates to asses that PK,t,E is "close" to p(- + E). The latter 
are based on the intuitive idea that t 2 ^fS^ I\ should be of order Kt 2 < G 2 /K when K is large 
and t < G/K. This can be proved using that for any K, T is stochastically dominated by a Pareto 
random variable Y distributed according to: 



r iipiic 
py(y) = < Q y2 



ify> \\p\U 

otherwise 



(A20) 



Such a random variable satisfies the following law of large numbers [43 

Y n = Yl + -- +y "4F 
n 



(A21) 
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where Y is distributed according to a stable law with exponent a — 1 and tail amplitude ||/?||oo- 
In particular, for all e > 0, a large enough and n large enough 

V(Y n > a) < (1 + e )ii__li___ (A22) 
a 

Henceforth, for t < G/K 

' (<" g r < - («r-i < A23 » 

for any a, 6 > 0, and if large enough (independent of E). This implies that vA- y X _ E ■ More 
generally, one can expand the distribution of T in powers of Kt 2 , which correspond to truncate the 
recursive equation (|7|) after a finite number of iterations. For the purpose of our work, it will be 
enough to deduce from (|A23[) that for all e > 0, a small enough and K large enough: 

Vt<G/K, Ve,E 6 1, inf p(e + e + E) - e < p K ,t,E( e ) < SU P P( e + e + E) + e (A24) 

e'S[-a,a] e'e[-a,a] 

Of course, the same holds with pK.t,E instead of px,t,E- In particular, we deduce from the above 
bounds that, for tK < G/K: 

lim sup |Afi:,t K ,,E(e) - p(e)\ = (A25) 

^ x ee[-m,m] 
\E\<1/K 
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